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Transition State Theory for dissipative systems without a dividing surface 
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Transition State Theory is a central cornerstone in reaction dynamics. Its key step is the identification 
of a dividing surface that is crossed only once by all reactive trajectories. This assumption is often badly 
violated, especially when the reactive system is coupled to an environment. The calculations made in this 
way then overestimate the reaction rate and the results depend critically on the choice of the dividing surface. 
In this Letter, we study the phase space of a stochastically driven system close to an energetic barrier in 
order to identify the geometric structure unambiguously determining the reactive trajectories, which is then 
incorporated in a simple rate formula for reactions in condensed phase that is both independent of the dividing 
surface and exact. 
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Transition State Theory (TST) provides the concep¬ 
tual framework for a large part of reaction rate the¬ 
ory. Originally developed to describe reactivity of small 
molecules^— , it was later extended to study a wide vari¬ 
ety of processes in very different fields that only have in 
common the existence of a transition from well-defined 
“reactant” to “product” states^ - —. Its great success 
comes from its simplicity, since it gives a straightforward 
answer to the two central problems in chemical dynam¬ 
ics: the identification of the reaction mechanism and a 
simple approximation to the reaction rate. 

The rate limiting step in many reactions is the crossing 
of an energetic barrier, the top of which forms a bottle¬ 
neck that the system must cross as the reaction takes 
place. If a dividing surface (DS) is placed close to this 
bottleneck, reaction rates can be readily computed from 
the steady-state flux through it. The TST approxima¬ 
tion is obtained under the assumption that the reactive 
classical trajectories cross the DS only once and never re¬ 
turn. Very often, for example when the system is strongly 
coupled to a noisy environment such as a liquid solvent, 
this no-recrossing assumption fails, and any conventional 
DS is crossed many times by a typical trajectory. As 
a result, TST calculations significantly overestimate re¬ 
action rates, and considerable effort has been devoted 
to the construction of a DS that minimizes recrossing^. 
The problem mainly derives from the fact that this spe¬ 
cial surface plays a double role: It defines (or separates) 
reactant and product regions, and it also serves to iden¬ 
tify the reactive trajectories. To the latter purpose the 
DS is ill suited. Despite this fundamental drawback, TST 
remains attractive since reactive trajectories are simply 
identified as those crossing the DS. This criterion only 
takes into account the instantaneous velocity at the DS, 
without requiring any time consuming trajectory simula¬ 
tion. 

Important advances in TST have recently been 
achieved within the approach of modern nonlinear dy¬ 


namics: A strictly recrossing free DS can be con¬ 
structed in the phase space of reactive systems with 
arbitrarily many degrees of freedom^iiy^-. These re¬ 
sults were generalized to systems interacting with an 
environment— r—. It has been shown that the desired 
exact TST can be constructed, in the harmonic limit, 
by using a moving DS only crossed once by the reac¬ 
tive trajectorie s 13 ’ 14 . Accurate results are still obtained 
for moderate anharmonicitieit and can be improved by 
normal form calculation s 18 : 19 . 

In this Letter, we make a step forward by presenting 
an explicit description of geometric phase space struc¬ 
tures in an anharmonic noisy system and an analytical 
scheme that relies on these structures for the compu¬ 
tation of exact TST reaction rates for arbitrary multi¬ 
dimensional potentials coupled to a noisy environment. 
The key point is the demonstration that reactive trajec¬ 
tories can be rigorously identified solely from their initial 
conditions, thus avoiding the choice of a (arbitrary) DS. 
This is done in terms of the stable manifold associated 
to a “noisy” Transition State (TS) trajectory jiggling in 
phase space. This geometric structure encodes the rele¬ 
vant information about the noise in the most economical 
manner and it can easily be incorporated into a rate cal¬ 
culation. Our method retains the fundamental simplicity 
of TST, also providing the conceptual tools to develop 
new computational algorithms. An application to the 
simple case of the one-dimensional quartic potential is 
presented as an illustration. 

The Langevin equation (LE) has been widely used 
to model the interaction of a reactive system with a 
surrounding heat bath22r— . Being a classical model, 
this description neglects quantum effects such as bar¬ 
rier tunnelling, which can be important in the case of 
light particles^, and the interaction with excited sur¬ 
faces through conical intersections— 

The dynamics of a unit mass particle defined by coor¬ 
dinate x moving in a one-dimensional potential U{x) is 
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given by 


X = ~U'{x) -J± + £ a (t). (1) 

Here, £ a (t) is the fluctuating force exerted by the bath, 
which is connected to the damping strength, 7 , by the 
fluctuation-dissipation theorem 

(Ut)Ut')) a = 2k B T'y8(t-t l ). (2) 



The model potential that we have chosen to study is 


U{x) 


1 

2 


L0lx 


2 
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although our derivation equally applies to any other one¬ 
dimensional case. This generality will be emphasized by 
using f(t), equal to — 04 a ; 3 in our case, to denote any 
anharmonic force. The extension to higher dimension is 
straightforward and will be presented elsewhere. 

For every fixed realization of the noise, the LE gives 
rise to a specific trajectory called TS trajector y 13 ' 14 . This 
orbit remains in the vicinity of the energetic barrier for all 
times, without ever descending into any of the potential 
wells. Other trajectories will be referred to the TS trajec¬ 
tory, which will be taken as a moving coordinate origin. 
Since the LE is a second order differential equation, its 
phase space is two-dimensional, with coordinates x and 
v = x. If we now introduce the new coordinates 


v — X s x 
A u — A s ’ 


s = 


—V + X U X 

A u — A s ’ 


with 


A 


U,S 


1 

2 


[~7± (7 2 + 4wb) 1/2 ], 
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relative values are defined by the time-dependent shift 
A u = u-u* J (t), As = s — s* (t), ( 6 ) 

where the coordinates of the TS-trajectory are 



The corresponding equations of motion are 

/( x* a + Au + As) 


Au = A U A u + 
As = A s As — 


Au A s 

/(a:* + Au + As) 
A u — A s ’ 


( 8 a) 

( 8 b) 


with a:* (t)=u£(f)+st (t) I and v i( t )=X u u t a (t)+X s s t a (t)]. 
Now, the geometric phase space structure in the vicin¬ 
ity of the barrier can be easily discussed. First, in the 
harmonic limit, i.e. f{x) = 0 , the equations of motion (151) 
are trivially solved, giving 

A u(t) = Au( 0 ) e Xnt , As(t) = As(0) e Xst . (9) 


FIG. 1. (color online) Schematic view of the phase 
space structure near the Transition State trajectory for the 
Langevin equation for the harmonic (a) and anharmonic (b) 
cases. The time-dependent invariant manifolds are attached 
to the TS trajectory, and move through phase space with it. 
In the harmonic limit they appear as (red) straight lines (a), 
but they get deformed by anharmonic couplings (b). Reactive 
and non-reactive trajectories are represented in black. 


Since A u > 0 and A s < 0, Au increases exponentially 
in time, whereas As shrinks accordingly. More impor¬ 
tantly, the lines Au = 0 and As = 0 are invariant under 
the dynamics of the system; being the unstable and sta¬ 
ble manifolds of the origin, respectively. As shown in 
Fig. HI a), these invariant manifolds separate trajectories 
with different qualitative behavior: Those above the sta¬ 
ble manifold (larger relative velocity) move to the prod¬ 
uct side in the distant future, while those below it move, 
on the other hand, towards the reactant side. 

In space fixed coordinates the invariant manifolds ap¬ 
pear attached to the TS trajectory, as shown by the 
dashed lines in Fig. [ljb); their instantaneous position 
depends on the realization of the noise. Accordingly, 
the manifolds move through phase space but they still 
separate trajectories with different asymptotic behaviors. 
The stable manifold intersects the x = 0 axis in a point 
with velocity V*. Trajectories with initial position x = 0 
and initial velocity v larger than this critical velocity, 
Vj, are reactive, while trajectories with initial velocities 
v < Vj are not. The (random) value Vj therefore en¬ 
codes the relevant information about the realization of 
the noise concisely. In other words, once the instanta¬ 
neous position of the stable manifold is known, any tra¬ 
jectory can unambiguously be classified from the values of 
its initial condition as reactive or non-reactive. Finally, 
the presence of anharmonicities (f(x) ^ 0 ) will distort 
the invariant manifolds, as indicated by the red lines with 
arrows in Fig. [lib). The main step of the theory to be 
developed here is the calculation of this deformation. 

This critical velocity can be calculated from the con¬ 
dition that the trajectory with a:(0) = 0 and u(0) = Vj 
is contained in the stable manifold of the TS trajectory. 
To find this trajectory, which will be called critical tra¬ 
jectory, we formally solve the equations of motion (JSJ by 


A u(t) = C u e Xut 


£ [A u , f(x* + Au + As); t) 
A u ~ A s 


( 10 a) 
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A s{t) = C s e Xst 


S[X Sl f(x$ + Am + As); t] 


A u — A 


S 


(10b) 


where C u and C s are two arbitrary constants, and the 
integral operator 


S T [g,g;t} 


/ OO 

g{r) exp [g(t — t)] dr : Reg > 0, 


+ / g( T ) exp [git - r)] dr 

Jo 


Reg < 0. 

( 11 ) 


has been introduced as a convenient shorthand notation. 
The subscript r indicating the integration variable in 
Eq. (fill) will be left out whenever this does not cause any 
ambiguity. The unknown constants C u and C s can be de¬ 
termined by noticing that the critical trajectory should 
approach the TS trajectory for large times. In particu¬ 
lar, it should remain bounded for t —> oo. This can only 
be satisfied if C u = 0. [This condition also ensures that 
the S functional in Eq. (llOal) is well defined.] With this 
choice, Eq. (IIOal) determines the initial condition Au(0), 
and C s = As(0) can then be found from the condition 
that x(0) = *1(0) + Au(0) + As(0) = 0. Finally, the 
coordinate transformation © yields 


V* = v(0) = (A u - As)'u(O) (12) 


In general, Eq. m represents only a formal solu¬ 
tion, since its right-hand side depends on the unknown 
functions Au and As. However, since it is known that 
the critical trajectory remains in the neighborhood of 
the barrier at all times, the anharmonic force will be 
small, and then its influence can be evaluated through 
perturbation theory, thus obtaining an expansion = 
Vq + C 4 V/ + C 4 V 2 + ■ ■ ■ in powers of the anharmonic 
coupling parameter C4. 

In the harmonic approximation the critical trajectory 
is given by 

A uo(t) = 0 and Asq(£) = —x :|: (0)e As *, (13) 


and for this case Eq. m yields 

V* = V* = (A u - A s )u*(0) 

that was already derived in Ref. 0. When the solu¬ 
tion (fl3l) is substituted into Eq. m, x = x$ + Au + As 
is replaced by 

x c*(t) = x i(t) ~ e Xst xi{ 0), (14) 


which is the harmonic approximation to the coordinate 
x(t) of the critical trajectory. Equation (II Oal) then gives 

A Ul (t) = -^— s[\ u j(xy,t], 

and therefore the leading-order velocity correction is 

V} = S[X u J(X);0] = -c 4 5[A u ,X3;0] (15) 



FIG. 2. Critical velocity for one realization of the noise as 
a function of the anharmonicity for cab = 1, 7 = 2.5, and 
ksT = 1: Numerical simulation results (red crosses), har¬ 
monic approximation (gray line), perturbative results to first- 
order (green straight line), and second-order (blue line). 


for the quartic potential ©. Given the initial condition 
Asi(0) = 1*4(0) — Aui(0), Asi(£) can be obtained from 
Eq. (IIObi) , and then substituted into Eq. (IIOal) to find 
Au 2 (t). In this way, the second-order velocity correction 


2 A„ - A s 


Sr 


X u,X 2 (r) (e AsT S'[A u , AT 3 ; 0] 


-S[X u ,X 3 -,t\ + S[X s ,X 3 ;t}) ;0 (16) 


is obtained. In Fig. [2] we present a comparison between 
numerical (exact) results for the critical velocity, Vj, for 
one realization of the noise and the approximate values 
computed using (fl5l) and m, as a function of the an¬ 
harmonic parameter, c 4 . 

In order to calculate the corresponding reaction rate, 
k , we choose the simplest DS, defined by x = 0, and use 
the basic flux-over-population rate formula (see e.g.— ), 
which states that k is proportional to the reactive flux 

j=(vy:t ( i? ) 

across the DS. This flux is to be averaged over different 
realizations of the noise, a, and also over a Boltzmann 
ensemble of initial velocities, v, for trajectories starting 
at the DS. Notice that only reactive trajectories should 
be included in the average. 

The non-recrossing assumption of conventional TST 
can be restated here by saying that the reactive trajec¬ 
tories are those crossing the DS with velocity v > 0. We 
call the rate constant obtained with this approximation 
k TST . Any effects beyond TST are customarily summa¬ 
rized^ into a transmission coefficient re = k/k TST . 

In terms of our stochastic invariant manifolds, reactive 
trajectories are characterized by v > Vj, as discussed 
before. Using this criterion, the Boltzmann average over 
velocities in Eq. 1171) can be evaluated, as it was in Ref. Util 
for the harmonic case. This gives the exact expression 


re = 



vi 2 


2k B T 


(18) 
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where only the average over the noise remains to be 
done. By substituting the perturbative expansion for 
the critical velocity into this expression and expanding 
the exponential, a perturbative series of rate corrections, 
k = Ko + C 4 K 1 + C 4 K 2 + ..., is obtained, where 

K o = ( E ) a . ( 19a ) 

= < I9b > 

with the abbreviated notation 


numerical scheme to achieve an efficient rate calculation 
in complex systems. 



E = exp 



= exp 


(A u -A s ) 2 4 2 (0) ~ 
2 k B T 


• ( 20 ) 


FIG. 3. Transmission factor as a function of the anharmonic- 
ity for the same values and legend of Fig. [2] 


Expressions similar to in can be obtained for the 
higher-order corrections. The leading order kq was eval¬ 
uated in Ref. GS It yields the well-known Kramers result 
for the transmission factor. With the result (fl5l) . the first 
correction term is given by 

ki = C4(A , U ~ As) ^ [A u ,(i?4(0^( T )) a ;0] . (21) 


To perform the remaining average and subsequently eval¬ 
uate the S functional, note that it* (i) and s* (i), and con¬ 
sequently X a (t), are Gaussian random processes, whose 
correlation functions were given in [TtL Details of this cal¬ 
culation are irrelevant for the purpose of this Letter and 
will be presented elsewhere. The final result is given by 


3 C 4 kBT f 1 — p 2 x 2 

«i = “7-— I 1 T7 —2 

4 V1 + M 


( 22 ) 


in terms of the dimensionless parameter p = A u /wb- This 
exp ression agrees with the corrections given in Refs. Ei- 
1271 . A comparison between perturbation theory and the 
results of a numerical simulation is shown in Fig. [3] along 
with the second-order perturbative correction, which can 
be obtained in a similar manner. The perturbative re¬ 
sults describe the rate correctly as long as the coupling 
is not too strong. For large negative values of C 4 the 
second-order correction loses its accuracy. By contrast, 
it is accurate for all positive C 4 shown in the figure. If C 4 
is increased further, the wells of the model potential @ 
become too shallow for a rate theory to be meaningful. 

In summary, we have demonstrated that the use of 
stochastic invariant structures allows a rigorous TST re¬ 
action rate calculation in general anharmonic and noisy 
systems. In our approach, the DS is only used to define 
reactant and product regions. To identify reactive trajec¬ 
tories we employ the stable manifold that is determined 
by the dynamics of the system itself. The resulting rate 
formula HID is not only remarkably compact, but also 
exact. In particular, the arbitrariness that is usually in¬ 
troduced by the choice of a DS is absent. Although our 
presentation is based on an analytical perturbation ex¬ 
pansion, the method can easily be incorporated into a 
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